On the absence of a spiral magnetic order in Li 2 Cu02 
with one-dimensional CuC>2 ribbon chains 
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On the basis of first principles density functional theory electronic structure calculations as well 
as classical spin analysis, we explored why the magnetic oxide Li2Cu02, consisting of CuCh ribbon 
chains made up of edge-sharing CuCU squares, does not exhibit a spiral-magnetic order. Our 
work shows that, due to the next-nearest-neighbor interchain interactions, the observed collinear 
magnetic structure becomes only slightly less stable than the spin-spiral ground state, and many 
states become nearly degenerate in energy with the observed collinear structure. This suggests 
that the collinear magnetic structure of Li2Cu02 is a consequence of order- by-disorder induced by 
next-nearest-neighbor interchain interactions. 

PACS numbers: 75.25.+Z, 75.10.Hk, 75.10.Pq, 71.70.Gm 



Copper oxides with C11O2 ribbon chains made up of 
edge-sharing CUO4 squares have one-dimensional chains 
of spin-i Cu 2+ ions, and exhibit unique physical proper- 
ties. LiCu^O^ and LiCuVO^— show ferroelectricity when 
their Cu02 ribbon chains undergo a spiral-magnetic or- 
der at low temperatures. For a chain of spin- 1 ions, a spin 
spiral structure is predicted when the nearest-neighbor 
(NN) ferromagnetic (FM) spin exchange J\ and the next- 
nearest-neighbor (NNN) antiferromagnetic (AFM) spin 
exchange J2 satisfy the condition IJ2/J1I > 0.25, while 
an FM structure is predicted if | J2/ J\\ < 0.25.— The cop- 
per oxide Li2Cu02 also consists of Cu02 ribbon chains, 
but has a different magnetic structure. A neutron powder 
diffraction study of Li 2 Cu02 at 1.5 K showed a collinear 
magnetic structure in which the spins of each Cu02 chain 
has an FM arrangement with Cu moments perpendicu- 
lar to the plane of the Cu02 ribbon and the arrange- 
ment between adjacent FM chains is AFM^ (hereafter 
this magnetic structure will be referred to as the AFM-I 
state). Thus, to explain this collinear magnetic structure, 
one might expect IJ2/J1I < 0.25 for the CuC>2 chains of 
Li2CuC>2. Indeed, Graaf et al. obtained IJ2/J1I = 0.15 
on the basis of first principles electronic structure calcu- 
lations using the embedded cluster model, 5 However, the 
Cu02 ribbon chains of Li 2 Cu02 are similar in structure 
to those of LiCu202 andLiCuVC-4, so that \ J2/J\\ > 0.25 
would have been expected. If \J2/Ji\ > 0.25, one needs 
to ask why a spiral magnetic order does not occur in 
Li2Cu02- In addition, more than two spin exchange in- 
teractions are necessary to describe the magnetic struc- 
ture of L12CUO2, and the nature and magnitude of these 
interactions are not unequivocal*^ 7 - Another puzzle con- 
cerning Li2Cu02 is that it undergoes a phase transition 
below ~2.4 K to a state believed to be a spin canted 
statei£&i£ So far, the origin and the nature of this phase 
transition remain unclear. 

The spiral magnetic order of LiCu202 and LiCuV04 
is a consequence of the spin frustration associated with 
the NN FM and NNN AFM interactions in their Cu0 2 
chains. A collinear magnetic order can occur as a conse- 
quence of order-by-disorder ) 1 1 i 12 which occurs typically 



in highly spin frustrated systems^ Provided that a spin 
spiral state is the ground state for the CuC>2 chains of 
Li2Cu02, one might speculate if the AFM-I state of 
Li2Cu02 is close in energy to the spin spiral state and if 
Li2Cu02 has a large number of nearly degenerate states 
around the AFM-I state. In the present work we explore 
these possibilities by studying the magnetic structure of 
Li2Cu02 on the basis of first principles density functional 
theory (DFT) electronic structure calculations and car- 
rying out classical spin analysis with the spin exchange 
parameters deduced from the DFT calculations. 

TABLE I: Relative energies (in meV/Cu) of the various mag- 
netic states with respect to the FM state obtained from 
GGA+U calculations with different U e ff values. 

Ueff (eV) 2 4 6 8 10 
E(AF1) -9.59 -5.80 -3.84 -2.38 -1.46 -0.96 
E(AF2) 0.69 4.89 5.76 5.82 4.70 3.67 
E(AF3) -15.34 -8.01 -4.29 -1.960 -0.71 -0.18 
E(AF4) -6.66 -1.00 1.23 2.23 2.22 1.88 
E(AF5) 1.25 5.00 5.77 5.80 4.65 3.65 



Our DFT electronic structure calculations employed 
the full-potential augmented plane wave plus local or- 
bital method as implemented in the WIEN2k code For 
the exchange-correlation energy functional, the general- 
ized gradient approximation (GGA)i^ was employed^ 
with BJ^^K max — 7.0. To properly describe the strong 
electron correlation in the 3d transition-metal oxide, the 
GGA plus on-site repulsion U method (GGA+U) was 
employed. 17 We also examined the energy of Li2Cu02 as 
a function of the magnetic order parameter q by employ- 
ing the non-collinear magnetism code, WIENncm^ 

Li 2 Cu02 has a body centered orthorhombic structure 
(space group Immm with a — 3.654 A, b = 2.860 A, and 
c = 9.377 A)^ where the CUO2 ribbon chains run along 
the b-direction (Fig.[lJ. As depicted in Fig.QJ,, there are 
five possible spin exchange interactions to consider; J\ 
and J2 are NN and NNN intrachain interactions, respec- 
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tively, J3 and J4 are NN and NNN interchain interactions 
along the c-direction, respectively, while J5 is the inter- 
chain interaction along the a-direction. To evaluate the 
interactions J\ — J5 , we calculate the relative energies of 
the six ordered collinear spin states shown in Fig. [2] in 
terms of GGA+U calculations. To see the dependence of 
these spin exchange interactions on the effective on-site 
repulsion U e ff — U — J, our GGA+U calculations were 
carried out with U e ff ranging from to 10 eV. (For 3d 
transition metals, U is generally less than 10 eV and the 
J value is usually 1 eV.) The relative energies of the six 
ordered spin states of Fig. [5] obtained from our GGA+U 
calculations are summarized in Table [IJ In terms of the 
exchange parameters J\ — J5, the energies of the six mag- 
netic states per Cu are written as 



E(FM) = (Jr + J 2 + 4J 3 +4J 4 + J 5 )/4 

E{AF1) = ( J l + J 2 - 4J 3 - 4J 4 + J B )/4 

E(AF2) = (- ,h + J 2 + J 5 )/4 

E(AF3) = (- J 2 + 2J 3 - 2J 4 + J 5 )/4 

E(AFA) = (-Ji+J 5 )/8 

E(AFb) = (-Ji + J 2 - J B )/8 



(1) 



Thus, by equating the energy differences of these states 
in terms of the spin exchange parameters with the cor- 
responding energy differences in terms of the GGA+U 
calculations, we obtain the values of J\ — J5 summarized 
in Table [Til where we employed the convention in which 
positive and negative numbers represent AFM and FM 
interactions, respectively. J5 is very weak in agreement 
with Mizuno et al.Z The NNN interchain interaction J4 is 
much stronger than the NN interchain interaction J3, and 
this finding does not support the assumption by Mizuno 
et al. that J3 and J4 are similar ii J4 is stronger than J3 
because the overlap between the magnetic orbitals, which 
depends on the overlap between the O 2p orbitals of the 
magnetic orbitals^ is much more favorable for the path 
J4 than for the path J3 (Fig. [3]). The NN intrachain in- 
teraction Ji is FM while the NNN intrachain interaction 
J 2 is AFM. These intrachain interactions are the same 
in nature to those reported by Graaf et a/.£ However, 
our study shows that | J 2 /Ji| > 0.25, for all U e fj values 
employed, and hence Li 2 Cu0 2 should have a spin-spiral 
ground state as far as isolated Cu0 2 ribbon chains are 
concerned. 



TABLE II: Calculated exchange parameters (in meV) de- 
duced from GGA+U calculations. 
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To see how the above prediction is affected by the 




FIG. 1: (Color online) (a) Crystal structure and five spin ex- 
change paths Ji — J5 of Li2Cu02. (b) Cu moments of the 
spin spiral ground state at q = (0, 0.20, 0) obtained from the 
GGA+U non-collinear calculation with U e ff = 6 eV. (c) De- 
tailed view of the Cu and O moments of a CUO2 ribbon chain 
in the spin spiral ground state shown in (b) . For the purpose 
of illustration, the O moments were increased by three times. 



(a) FM 



(b) AFM! 



(c) AFM2 



4b 




(d) AFM3 



(e) AFM4 



(t) AFM5 



4b 





FIG. 2: (Color online) Schematic representations of the six 
ordered spin arrangements of Li2Cu02 employed for GGA+U 
calculations to extract the five spin exchange parameters Ji — 
J5. The filled and empty circles the up-spin and down-spin 
Cu sites, respectively. 
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NNN interchain interaction J4, we carried out a clas- 
sical spin analysis based on the Freiser metho d 20 ! 21 using 
the three dominant exchange parameters Ji, J2, and J4. 
The spin interaction energy of an ordered spin state with 
q = (2irq x /a, 2nq y /b, 2-Kq z /c) can be written as 

E(q) = J 4 {cos[27r( & /2 + 3 fe /2 + <fe/2)] 
+cos[2ir(-q x /2 + 3q y /2 + q z /2)} 
+cos[2w(q x /2-3q y /2 + q z /2)} (2) 
+cos[27r(( ?:c /2 + 3q y /2 - q z /2)}} 
+cos(27T(j a ) Ji + cos^^g,,) J 2 . 

For simplicity of our discussion, we will represent q by 
(q x , q y , q z ). This E(q) vs. q relation has minima along 
the (0, q y ,0) direction. The E(0, q y , 0) vs. (0, q y , 0) 
curves calculated with the spin exchange parameters de- 
rived from the GGA+U calculations for U e f / = 6 eV are 
presented in Fig. 21 The solid curve, obtained only with 
the intrachain interactions J\ and J2, shows two minima 
(at q y = 0.18 and q y = 0.82) of equal energy. The FM 
state {q y = 0.00) and the AFM-I state (q y = 1.00) are 
identical in energy, and are less stable than the two spin- 
spiral states (q y = 0.18 and q y — 0.82). These are the 
expected results in the absence of the interchain inter- 
action because IJ2/J1I > 0.25. The dashed curve, ob- 
tained with the intrachain interactions J\ and J2 as well 
as the interchain interaction J4, also shows two minima 
at q y = 0.21 and q y = 0.90. Note that the interchain in- 
teraction J4 raises the energy of the FM state while low- 
ering that of the AFM-I state. As a result, the E(0, q y , 
0) vs. (0, q y , 0) curve around q y = 0.21 becomes sharper 
while that around q y — 0.90 — 1.00 is nearly flat. Both 
spin-spiral states are only slightly more stable than the 
collinear AFM-I state. Our calculations using the spin 
exchange parameters obtained with U e ff < 6 eV show 
that the energy around q y — 0.21 becomes lower than 
that around q y = 0.90, and both states have lower ener- 
gies than the collinear AFM-I state (q y — 1.00). In terms 
of the parameters obtained for U e ff > 6 eV, however, the 
collinear AFM-I state becomes the ground state. 

Now we evaluate E(0, q y , 0) vs. (0, q y , 0) relations on 
the basis of non-collinear GGA+U electronic structure 
calculations using the WIENncm code*^ In this method, 
the incommensurate spiral magnetic order is simulated 
without resorting to the supercell technique by using the 
generalized Bloch theorem^ The E(0, q y , 0) vs. (0, q y , 
0) relation calculated for the representative C/ e // (i.e., 6 
eV), presented in Fig.[l]as a dotted line, is quite similar to 
that found from the classical spin analysis. An important 
difference is that the non-collinear GGA+U calculations 
predict the spin-spiral state at q y = 0.20 to be slightly 
more stable that that at q y — 0.95. The spin arrange- 
ment of the spin-spiral state at q y — 0.20 is illustrated in 
Fig. [TJd and QJ;. In this state of zero total spin moment, 
the non-collinearity of the spin arrangement occurs not 
only between Cu spins but also between the O and Cu 
spins. Our calculations show substantial moments on the 
O sites, as found in the previous studies^ ' 10 ! 23 From our 



calculation with U e ff = 6.0 eV, the oxygen spin moment 
is 0.11 hb, which agrees with the LDA+U result^ and 
the experimental value (between 0.10 and 0.12 hb)*^ 
Our non-collinear GGA+U electronic structure calcula- 
tions with U e ff > 6 eV or with U e ff < 6 eV still show 
that the ground state is a spin-spiral state. Thus, with 
any reasonable U value, we predict a spin-spiral ground 
state for Li 2 Cu02. 




(a) J 3 (b) J 4 

FIG. 3: Arrangements of the C11O4 squares and their mag- 
netic orbitals associated with (a) the NN interchain interac- 
tion J3 and (b) the NNN interchain interaction J4. The two 
adjacent CUO2 ribbon chains differ in their a-axis heights by 
a/2. The CUO4 squares with different a-axis heights are indi- 
cated bv thick and thin lines. 
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FIG. 4: (Color online) E(0, q y , 0) vs. (0, q y , 0) relations cal- 
culated for Li2Cu02. The solid and dashed lines are based 
on the classical spin analysis (solid line: only with the intra- 
chain interactions Ji and J2, dashed line: with the intrachain 
interactions Ji and J2 as well as the interchain interaction 
J4). The dotted line is based on non-collinear GGA+U cal- 
culations with U e ff = 6 eV, where the circles represent the 
calculated points. 

From our non-collinear GGA+U calculations, the en- 
ergy difference between the spin-spiral state at q = 
(0,0.20,0) and AFM-I state at q = (0,1.00,0) is very 
small (Fig. 2]). In the case of U e f / = 6 eV, the difference 
is 1 meV/Cu and deceases with increasing U e ff. From 
the classical spin analysis shown in Fig. |H this energy 
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difference is even smailer. As aiready pointed out, the 
E(0, q y , 0) vs. (0, q y , 0) curve is sharp around q y = 0.20 
but nearly flat around q y = 0.90 — 1.00. As a conse- 
quence, the states in the region of q y = 0.90 — 1.00 are 
nearly degenerate, and are only slightly less stable than 
the spin-spiral ground state at q y ~ 0.20, namely, the 
density of states is much higher in the region of the AFM- 
I state than around the spin-spiral ground state. The lat- 
ter provides a natural explanation for why the CuC>2 rib- 
bon chains of Li 2 Cu02 do not exhibit a spiral-magnetic 
order despite that the CuC>2 chains are very similar in 
structure to those found in LiCu2C>2 and LiCuV04, and 
Li2Cu02 has a spin-spiral ground state. In short, the 
AFM-I structure (q y — 1.00) is a collinear order arising 



from the occupation of many nearly degenerate states 
around q y = 0.90 — 1.10, and hence is an example of 
order-by-disorder . 1 1 i 12 The phase transition below 2.4 K, 
believed to be a transition to a spin canted state, might 
arise from an increased population of the spin-spiral state 
at q = (0,0.20,0). What distinguishes Li 2 Cu02 from 
LiCu2C>2 and LiCuV04 is the NNN interchain interac- 
tion J4, which lowers the energy of the states around the 
AFM-I state and makes them nearly degenerate. 
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